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ABSTRACT 

A three-dimensional methematical model to predict velocities 
and temperature aneimolies is presented. The effects of 
bottom topography, wind, surface heat transfer and currents 
are included. A co-ordinate transformation which maps 
variable depth basin to a constant one is used. The rigid- lid 
assumption is invoked o A Poissons equation is derived as the 
predictive equation for surface pressure. Explicit schemes 
are used in time to integrate the primitive conservation 
equations and an iterative scheme is used to solve the 
Poisson's equation. The model was applied to Florida Power 
and Light (FPL) company's fossil fuel power plant at 
Cutler Ridge site for April 15,1975, outfall conditons. 

Infra red (IR) data in the morning (0912 EST) on April 15, 
1975 was used as initial conditions for temperature. The 
model vas run for 2 hours 38 minutes and compared with IR 
data at (1150 EST) on the same day. The predicted results 
compare favorably with IR data. 
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INTRODUCTION 

A major source of thermal pollutxon in rivers, lakes and 
estuaries is the surface discharge of heated water from 
condensers of fossil and nuclear power plants o In order to 
design intake and outfall structures to avoid recirculation, 
sufficient knowledge of flow and temperature distribution 
in the receiving water near the discharge is necessary. 

Most of the previous work for the near field problem is of 
integral type. Some of the important work in integral models 
are that of Carter and Reiger (1974), Hoopes, Zeller and 
Rohlich (1971), Stolzenbach and Harleman (1973), Hirst (1972) 
and Stefan and Vaidyaraman (1971). Dunn, Policastro and 
Paddock (1975) discuss most of the integral models and their 
drawbacks. They also discuss mamerical models, their 
advantages and disadvantages- Some of the numerical models 
for heated discharges are those of Brady and Geyer (1972), 
and Waldrop and Farmer (1974). All the above mentioned models 
both integral as well as numerical suffer some common dis- 
advantages. These models also do not have adequate verifica- 
tion. The handbook prepared by Shirazi and Davis (1974) for 
surface ischarges using formulations by Prych (1972) is the 
most readily usable predictive tool available at present. 

To date, there is no three-dimensional model which incor- 
porates the effects of bottom topography, wind, surface 
heat transfer and current which has been extensively veri- 
fied. These factors play an important role in the dispersion 
of heat. The present approach tries to incorporate all these 
relevant mechanisms using the divergence of Navier Stokes 
model developed by Sergupta and LicK (1974). 

For calibration and model verification, FPL’s power plant at 
Cutler Ridge site in Florida is considered. The objective 
of the present investigation was to predict velocity and 
temperature anamolies in the near field (i.e., the region 
which IS most affected by the heated discharge) . The 
region of interest is approximated to 525m along the axis 
of the jet and 425m across the axis of the jet. Figs. 1 
and 2 show the basin and grid system. 
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BASIC E':2UATI0NS 

The equations which describe the behaviour of fluid flow are 
those expressing the conservation of mass, momentum and 
energy. In order to keep the generalized nature of the model 
all significant terms in the conservation equations are re- 
tained. This includes the effects of buoyancy, inertia, 
coriolis, density and turbulent mixing. Wind shear, surface 
heat transfer and cross flow are also considered. The details 
of the derivation of the equations are presented in Sengupta 
and Lick (1974) . 


Assumptions and Approximations 

The riqid-lid approximation is made. The riaid-lid approximation 
IS that a rigid wall is placed on the top which allows hori- 
zontal velocities, but forces vertical velocities at the sur- 
face to zero. This approximation was used by Bryan (1969) 
Sengupta and Lick (1974) and several others. 


This approximation is permissible as it eliminates surface 
gravity waves and are therein associated instabilities. This 
approximation allows use of larger time steps thus saving 
computer time. 

Since the velocities are small, the vertical momentum equation 
IS replaced by the hydrostatic approximiation. 

Turbulence is modelled using eddy transport coefficients. 
Constant, but different values are used for horizontal and 
vertical directions. The eddy transport coefficients are 
found by trial and error method till the correct values are 
obtained. 

In order to overcome programming difficulties the equations 
are stretched vertically using the transformation, 

V =Z A (x,y) 

Now the equations will be transformed from (x,y,a) coordinate 
system to a more convenient coordinate system (x,y,v)or for 
convenience (cif,p,Y)- The advantage of this transformation is 
that the same number of grid points can be used in the vertical 
direction at shallow as well as deeper parts. Now the equa- 
tions after using the above transformation will take the form 
snown re low: 
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The last term in equation (7) is the Hirt and Harlow (1967) 
correction term. This term is not set equal to zero even 
though the rigid- lid conditions imply W (S = 0) =0. This 
IS in anticipation of the errors which will arise when the 
Poisson equation is solved by iterative technique. If 
these errors are not corrected, the continuity equation 
will not be satisfied leading to acc\imilation or loss of 
fluid from the system. 


Equations 1 through 7 form the set of equations which 
are to be solved with appropriate boundary conditions. 
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I 

Bo\indarv Conditions 

; The "boundary conditxons used are shown in Fig. 3- The condi- 

tions on solid walls and bottom are no slip and no normal 
velocity for all time for the momentum equations, except 
W / 0 at side wall due to hydrostatic approximation. The 
temperature boundary condition at solid walls and bottom 
IS handled by assuming the walls and bottom adiabatic. The 
boundary condition for velocity and temperature at open 
boundaries is ^f^ = 0, where ^_f is velocity or temperature 

gradient normal to the boundary. At the air-water interface 
the boundary condition on the momentum equations are obtained 
, in the form of wind stress. The temperature boundary condi- 

tion at the air water interface due to Edinger & Geyer (1971)is 

Heat flux = AQ = K (T„- T„) 

S -bi w 

where Ks is surface heat transfer coefficient and T_ and T_ 
are equilibrium and surface temperatures. The boundary condi- 
tions for the Poisson equation are Neuman type along solid 
’ boundary and along open boundaries the pressure is assumed 

constant. 

The boundary conditions in summary for the stretched system 
of aquations are: 

Along the shore: 
u = 0 

V = 0 

n / 0 

At the bottom of the basin (y=1) : 

n = 0 

u = 0 

V = 0 
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Alonq frca boundaries: 
At (of = 0,P,y) 

^1=1, K =^I=2,K 
V = 0 
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m m 
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n = 0 (Rigid-lid) 
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ARRANGEMENT OF VARIABLES AND GRID SYSTEM 


The grid system and arrangement of variables are shown in Pig. 
4. u and v are located on integral nodes (I,J,K) on the hori- 
zontal (x,y) plane. P,T,W,p are located at half grid points 
{1 + hi J hi k) . This arrangement is repeated in Z direc- 
tion. This kind of arrangement allows better meshing of the 
variables without calculating all variables at all nodes of 
the staggered mesh system. Constant grid spacing is used in 
a or a direction. Constant grid spacing is also used in x 
and y (or a and p) directions. 


METHOD OF SOLUTION 

The method of solution is explicit using forward time central 
space (FTCS) schemes except in the diffusion term. The 
diffusion terms are differenced by the Dufort-Frankel scheme. 
The initial values of variables are specified and then the 
equations are integrated marching in time. The only itera- 
tive procedure involves the solution of the equation for 

surface pressure. 

Since it is not possible to make a strict stability analysis 
for the system of equation under consideration, one-dimen- 
sional Burgers equation is used for stability analysis. 

The stability ciriterion for the equations under consideration 
is therefore: 
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The equations are solved on a UNIVAC/1106 computer at the 
University of Miami. 

RESULTS 

The region of influence of the Cutler Ridge plume has been 
approximated to a rectangular domain as shown in Fig. 1. 

It extends 425 meters in lateral extent and 525 meters in 
longitudinal extent. The discharge is taken as 25 meters 
wide. The numerical grid system is shown in Fig. 2. It 
has 18 and 22 nodes across and along the axis of the jet 
respectively. There are 5 nodes in the vertical direction. 

In order to understand the physical processes involved and 
to investigate the numerical behaviour of the model several 
simplified cases were executed before the final calibra- 
tion and verification run for April 15,1975. 


Case 1 


This case was designed to study the effects of bottom topo- 
graphy. A section of bottom topography can be seen in 
Fig. 5. The topography is symmetrical about the discharge 
centerline. The canal is 6 meters deep gradually merging 
into the bay which is 1.2 meters deep. The steady state 
surface velocities are symmetrical as expected, since the 
effects of wind and current are neglected and the bottom 
topography is symmetrical about the centerline. These 
results are shown in Fig. 6. There is a reduction in velo- 
cities in the bottom layers owing to bottom shear. Fig. 5 
shows centerline surface velocity with distance from dis- 
charge origin at different times. The steady state curve 
(155 min) shows that the velocity decays to about 75 meters 
and then gradually increases to a value slightly higher 
than the initial discharge velocity at a distance of 175 
meters. The reason for this is that the depth is initially 
constant along the jet and the velocity decays as expected. 
However, after 75 meters the depth decreases linearly to 
1.2 meters in the next 75 meters thus causing a constrict- 
ing effect which increases the velocity offsetting the 
decrease caused by lateral mixing. Here again it can be 
seen that the time to reach steady state increases with dis- 
tance from the jet origin. 


- 12 - 



VI-B-97 


Case 2 


The effects of density variations is studied in this case. 
The basin was taken as constant depth 1.2 meters and effects 
of wind^and current neglected. The discharge^ temperature 
IS 35.9 C and the ambient temperature is 28.0 C. Fig. 7 
shows tho surface isotherms. The .inertial core of the dis- 
charge can be seen as an elongated wedge of warm water at 
the mouth of the discharge. The isotherms are symmetrical 
as expected, since the effects of wind and current are neg- 
lected. The centerline temperature decay of the plume is 
shown in Fig. 8. The behaviour is very similar to the sur- 
face velocity decay. This is because the dispersal of heat 
and momentum essentially have the same mechanisms. The 
temperature decay near the origin is therefore followed by 
a region of rapid decay caused by entrainment. Finally as 
the difference in temperature is reduced the temperature 
asymptotically approaches the ambient value. 


Case 3 


This case includes the effects of strong wind along with 
density variations on a discharged plume into a constant 
depth basin. The wind was blowing at 6.71 m/sec (15 mph) 
160°N. Fig. 9 shov;s surface velocities. The wind driven 
current is very dominant for most of the flow. The effect 
of wind slowly reduced in the bottom layers, finally being 
dominated by bottom shear near the floor of the basin. 

Fig. 10 shows the skewed plume. The plume is turned towards 
the west. The skewing increases with time until a steady 
state is reached. 


Case 4 


The data base for this case is obtained from the field 
experiments conducted at the Cutler Ridge site on April 15, 
1975. The initial temperature conditions are taken from 
the morning IR data (0911-0912 EST) and ground truth data 
on April 15,1975. The computations were continued for 
2 hrs. 38 min. Fig. 11 shows isotherms predicted by the 
model along with IR data (11:45-11:55 EST) on April 15,1975. 
As can be seen there is reasonably good agreement between 
IR data and model predicted results. Fig. 12 shows tempera- 
ture decay along J at 1=11 (close to centerline) predicted 
by the model along with IR data. There is a good agreement 
between IR data and model results. The isotherms along I 
and J sections are shown in Figs. 13 and 14 respectivelv. 
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TABLE 1 


CASE NO. 

DEPTH 

DENSITY 

WIND STRESS 
(WIND SPEED) 

CURRENT 

REMARKS 

1 

Variable 

Constant 

None 

None 

Steady State 
(TTOT =155 min) 
At = 30 sec. 

2 

Constant 
(1.2 m) 

Variable 

None 

None 

TTOT = 1/2 hr. 
At = 15 sec. 

3 

Constant 
(1.2 m) 

Variable 

15 mph 
6.71 m/sec 

None 

TTOT = 2 1/2 hrs. 
At = 15 sec 

4 

Constant 
(1.2 m) 

Variable 

15 mph 
6.71 m/sec 

3 cm/sec 
south 

TTOT = 2 hrs. 38 min 
At = 10 sec. 


The following values are used in 
Discharge Vel ; 20 cm/sec 


Discharge Temp 
Initial Temp 
Air Temp 


35,9°C 

28.0°C 

29.5°C 


the above cases : 

K = 100 BTU/°F-Ft^-day 

O 


Case 4 is ran with IR data initial conditions 
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Stratification near discharge with isothermal conditions 
away from the plume can be seen in Fig. 13. 


SUMMARY 

Simple cases have been investigated to obtain an understanding 
of the basic mechanisms of the behaviour of a buoyant surface 
jet. Effects of wind, topography and current have been 
documented o Comparison with data bases obtained for field 
measurements, (both in-situ and airborne IR) is good. The 
model can be used either for preparation of environmental 
impact satements or research in buoyant surface jets. 
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HORIZCW'J'AL PlAUi: 



Arrangement repeated at each horizontal level 



Fig. 4 Arrangement of Staggered 
Grid and Variables 
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Fig« 5 Center Line Velocity Decay (Rigid-Lid Model) (Case 1) 
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Fig. 6 Surface Velocity Distribution 
at Steady State For Cutler 
Ridge Site (Rigid-lid Model) 
(case 1) 


20 cm/r.oc 
Var j able 
None 

Constant 


50 m 3 no m 


- 22 - 



VI-B-107 


DlscharRO Velocity 

Discharge Tenperaturo 

Depth 

Dcnsl ty 

ft' 1 nd 

“Talr 

Tinltial 

^total 


20 cm/soc 
35, 9 C 
1.2 m 
Variable 
None 
20. G°C 
28.0°C 


30 nln. 


CL 


o o o o 

o o o o 
lO o o 

• • • • 

00 a\ a cn 

cj M wo 


o 

td 

u 


3 « 

HU Ij 


■ ‘ V I ; ; / « 

« I ‘ “ 
I i • j 
‘ / / / ;; 

I I / / 

' ' ; ' 

1 1 , j 

\ ‘ \ '' ' '■ 

\ W\ \ 


oo o o o o 

o o o o o o 

u> O O lO O lO 

ift d- CO o <n CO 

CO^ CO CJ CVI CJ 


II 1 

' 1 • • ■' 
I ' 


li 


\ V \ \ V 
\ ^ \ \ 


, 1 . t 

I I ' » ' 

! !’■; 

f I f > 

/ ' f ' 

I /' / 




\ 

N N 
\ " 
> 

N. 


' ' / 

■'y'/ 




— 


0 

L 


25m 

J 


Fig« 7 Surface Isotherms for Cutler 
Ridge Site (Rigid-Lid) 

(Case 2) 
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Temperaturo ( °C ) 
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Di&char'tc Velocity 

Discliarfrc Temperature 

Depth 

Density 

Wind 
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,air 
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^lotal 


20 cm /see 
35, 9°C 
1,2 m 
Variable 
None 
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28,0°C 

30 min. 



Fig« 8 Centerline Temperature Decay For 

Cutler Ridge Site (Rigid-Lid Model) 
(Case 2) 
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Fig« 9 Surface Velocity Distribution 2h 
Hovurs After Discharge For Cutler 
Ridge Site (Rigid-Lid Model) 

(Case 3) 
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Fig« 10 Surface Isotherms for 
Cutler Ridge Site 
(Rigid-Lid) (Case 3) 
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Fig., 12 Comparison of Tenp Decay at I = 10 Along J 

(close to centerline) for April 15,1975 (11:55) 
(case 4) 
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Fig. 13 Vertical Section Isotherms Along Canal 
Center Line for April 15,1975 at Cutler 
Ridge Site (Rigid-Lid) (Case 4) 


- 29 - 




Depth (meters) Depth (meters) Depth (meter 


VI-B-114 



Fig o 14 Vertical- Section Isotiierms at Different 
J Sections for April 15,1975 at 
Cutler Ridge Site (Rigid-Lid) (Case 4) 
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